Dynamics of a liquid dielectric attracted by a cylindrical capacitor 
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The dynamics of a liquid dielectric attracted by a vertical cylindrical capacitor is studied. 
Contrary to what might be expected from the standard calculation of the force exerted by 
the capacitor, the motion of the dielectric is different depending on whether the charge or the 
voltage of the capacitor is held constant. The problem turns out to be an unconventional 
example of dynamics of a system with variable mass, whose velocity can, in certain cir- 
cumstances, suffer abrupt changes. Under the hypothesis that the voltage remains constant 
the motion is described in qualitative and quantitative details, and a very brief qualitative 
discussion is made of the constant charge case. 
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I. INTRODUCTION 

Electric forces on insulating bodies are notoriously difficult to calculate [1]. Only in very special 
instances can the calculation be performed by elementary means. The simplest case, treated in 
every undergraduate textbook on electricity and magnetism, is that of a dielectric slab attracted 
by a parallel-plate capacitor. Nearly as simple is the determination of the force of attraction on 
a liquid dielectric by a vertical cylindrical capacitor. It is remarkable that, under appropriate 
conditions, in both cases conservation of energy alone leads to a simple expression for the force 
with no need to know the complicated fringing field responsible for the attraction. Conservation of 
energy leads to the same expression for the force independently of whether the voltage or the charge 
of the capacitor is assumed to remain constant during a hypothetical infinitesimal displacement 
of the dielectric. It is reassuring that a detailed treatment based on a direct consideration of the 
fringing field leads to the same result Qj. Since the dynamics of the dielectric under such a force 
is not discussed in the textbooks, one may be left with the false impression that the motion is 
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unaffected by the choice of constant voltage or constant charge. 

In this paper we study the dynamics of a liquid dielectric pulled into a vertical cylindrical 
capacitor by electrostatic forces. The motion is quantitatively different, though qualitatively the 
same, depending on whether the voltage or the charge of the capacitor is held constant. It will be 
seen that this is a very peculiar variable-mass system, since in certain circumstances the velocity 
can change abruptly during a very short time interval right after the beginning of the motion. 

The paper is organized as follows. In Section 2 we review the elementary calculation of the 
force exerted by a vertical cylindrical capacitor on a liquid dielectric. In Section 3 the dynamics is 
investigated assuming the voltage applied to the capacitor stays constant. The motion of the liquid 
near the bottom edge is described qualitatively only, since our expression for the electrostatic force 
is not valid near the ends of the capacitor. For the motion when the surface of the dielectric is 
well inside the capacitor quantitative results are obtained. In Section 3 we briefly discuss the case 
in which the charge of the capacitor is held fixed, and a qualitative description of the motion is 
given. Section 4 is devoted to final comments and conclusions. 

II. FORCE OF A CYLINDRICAL CAPACITOR ON A LIQUID DIELECTRIC 

Let us consider a coaxial cylindrical capacitor of length L whose inner and outer radii are a 
and b, respectively, with L>t. Suppose the capacitor is in the vertical position and is partially 
immersed in a liquid dielectric contained in a tank, in such a way that the surface of the liquid 
that has entered the capacitor is far from its ends. If a potential difference is applied between the 
capacitor plates, the dielectric will be pulled into the capacitor. The calculation of the force of 
attraction is a textbook problem [3] which we proceed to work out for subsequent use. 

Let Lq be the length of the capacitor that is immersed in the liquid before the voltage is turned 
on. After the voltage is turned on the liquid inside the capacitor rises to the height Lq + x , as 
depicted in Fig 1. Let Ai and A2 be the linear charge densities on the parts of the inner plate 
that are in empty space and in contact with the dielectric, respectively. The corresponding electric 
fields in vacuum and inside the dielectric are 



Ei = - s , E 2 = s , 1 

Ineos lues 

where s is the radial unit vector and e is the permittivity of the dielectric. The voltage V applied 
to the capacitor is given by the radial line integral either of Ei or E 2 : 
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FIG. 1: A cylindrical capacitor partially immersed in a tank containing a liquid dielectric. The liquid inside 
the capacitor rises when a potential difference is applied between the plates. 
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V= I Ei -sds = — — In- = / B 2 -sds = — In-. (2) 
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From this the capacitance can be readily found. Indeed, the total charge of the capacitor is 



Q = X 1 [L- (L + x)} + A 2 (L +x) = Jj^{ eo [L - (L + x)] + e(L + x)} , (3) 

whence 



c = RW |i+Xe(Lo+:l01 ' (4) 

where Xe is the electric susceptibility of the dielectric. 

The force on the dielectric can be determined by equating the work done by the capacitor 
during a displacement dx of the dielectric to the decrease of the electrostatic energy stored in the 
capacitor. Thus, assuming the charge Q of the capacitor remains constant as the displacement 
takes place, 



dE d (Q 2 \ lQ 2 dC neoXeV 2 

(5) 



dx dx\2CJ 2C 2 dx ln(6/a) ' 



4 



where Eq. Q has been used. Since this force is positive, the dielectric is pulled into the capacitor. 
As is well known, the same result is obtained for the force if, instead of the charge Q, the voltage V 
of the capacitor is assumed to remain constant during the displacement. But in the latter case the 
battery responsible for maintaining the voltage fixed does work as the dielectric is displaced, and 
this work must be taken into account in the energy balance Q]. Since the expression for the force 
is the same in both cases, it may seem obvious that the motion of the liquid dielectric does not 
depend on whether the charge or the voltage is held constant. This is not the case, however, for a 
reason that may be easily overlooked. If the voltage is held constant then the force ([5]) is constant 
during the entire motion. However, if the charge Q remains constant then, using V = Q/C and 
Eq. (@|), the force ([5]) becomes 

Q 2 TrepXe _ Q 2 Xe Mb/ a) 

C 2 ln(b/a) 4vr6o[L + Xe (L + x)] 2 ' U 
This is a variable force that depends on the position of the dielectric because the capacitance 
increases and the voltage decreases as the dielectric penetrates the capacitor. Therefore, the 
dynamics of the dielectric is different depending on whether the charge or the voltage of the 
capacitor is constant throughout the motion. 



III. EQUATION OF MOTION OF THE LIQUID DIELECTRIC 

The equation of motion of the liquid dielectric is 



— (mx)=F + pgL Tr(b 2 -a 2 )-mg (7) 

where m is the mass of the dielectric that is already inside the capacitor at time t. The presence 
of the term proportional to L$ on the right-hand side of (|7|) is justified as follows. The pressure 
at the bottom of the liquid column that is inside the capacitor is higher than the pressure at the 
top by pgLo, where p is the density of the dielectric. This gives rise to a net upward pressure force 
that equals pgL§ir(b 2 — a 2 ) and must be added to the electrostatic force F. We have assumed that 
the area of the tank is much larger than ir(b 2 — a 2 ), so that the variation of the height of the liquid 
in the tank can be safely ignored. Other hypotheses will be mentioned in due time. 
The mass of liquid inside the capacitor is 



m = 7r(b 2 — a 2 )(Lo + x)p 



(8) 



5 



Insertion of this expression for the mass into the equation of motion ([7]) leads to 



Jt [iL ° + x)±] = ^-a^p - 9X - (9) 



IV. MOTION IN THE CONSTANT VOLTAGE CASE 



If V is constant, the use of allows us to write Eq.® as 



yy + y + gy = gi* + p{p2 _ ^ ■ (10) 



where we have introduced the new variable y defined by 



y = L + x. (11) 

The task of solving equation (jlOp is made easier by writing it in dimensionless form. It is easy 
to check that the positive constants ao and to defined by 

°° " pg(V-a*)]n(bM ~ ° ( ' 

have dimensions of length and time, respectively. In terms of the dimensionless variables £ and r, 
and dimensionless parameter a, defined by 



y L + x t L 

£ = — = , t = — , a = — , (13) 

ao ao to «o 

the differential equation ([TP]) reduces to 



' dr 2 V dr 

This equation has the equilibrium solution £ e(? = 1 + a, which is tantamount to x eq = oq. Thus, ao 
is the height at which the weight of the dielectric that has risen above the external level is exactly 
counterbalanced by the electric force that pulls it up. 

Equation (|14h can be solved by standard techniques. First of all we make a change of indepen- 
dent variable: 
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= d£ d?(j = dp = dpd£ = dp 
P ~ dr dr 2 ~ dr~ d£ dr ~ P d£ ' 1 j 



This reduces Eq.(|M]) to 



e^+P 2 + e-i-« = o, (i6) 



or, in terms of differential forms, 



tpdp+(p 2 + £-l-a)dZ = 0. (17) 

We seek an integrating factor iV(£) that turns the left-hand side of this equation into an exact 
differential: 



N(0£pdp + n(0(p 2 + £ - 1 - «K = • (is) 

Such a function G(£,p) exists if and only if 



Q^MQfr] = g^MOip 2 + H-l-a)]^Z— = N. (19) 
Thus, iV(£) = £ is an integrating factor and we have 



dG = fpdp + £(p 2 + £ - 1 - a)d£ . (20) 

The function G can be readily found by integrating dG along the broken path composed of the two 
straight line segments (0,0) — > (£,0) — ► (£,£>), and one immediately gets 



G(e,p) = ^y+f - l -^e. (2i) 

Therefore, the general solution of (|16p is 



^V + f -^£ 2 = C, (22) 



where C is an integration constant. 
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A. Close to the edge: a qualitative description 



Equation ([5]) is adequate for the force only when the surface of the liquid is far from the edges 
of the capacitor. Therefore, it is not legitimate to put Lq = in equation (llOp . Nevertheless, 
something remarkable happens when one considers Lq = 0. Since the attraction is caused by the 
fringing field, it is clear that the force on the liquid dielectric is not zero when its surface just 
touches the lower edge of the capacitor. Suppose that at t = the liquid is at rest, the voltage is 
turned on and the upward motion begins with £ = 0. Then the "first" ascending layer of the liquid 
has zero mass but is pulled into the capacitor with a finite force. Thus, ideally, it is subjected 
to an infinite acceleration, which brings about a finite jump of the velocity. Sure enough there is 
a transient process during which the electric field builds up and the velocity changes continuosly. 
However, since the time scale of this transient process is certainly much shorter than the time 
scale of the motion, it seems reasonable to expect that it can be satisfactorily modelled by a 
discontinuous jump of the velocity. Furthermore, the value attained by the velocity right after the 
transient, when so to speak the electrostatic force takes over, most likely cannot be freely chosen. 
This is a very peculiar behavior of this unconventional variable-mass system, whose mass can drop 
to zero at certain instants in the course of its motion. 

It turns out that this unique behavior is nicely illustrated, at least in some of its main qualitative 
aspects, by our description. In order to check this statement let us try to find out what the equation 
of motion has to tell us about the dynamics of the liquid near the edge of the capacitor, although 
the quantitative results are admittedly unreliable. 

With a = because we are taking Lq = 0, the initial condition £(0) =0 demands p(0) = ±1, 
as plainly seen from either (|14p or (|16p . Accordingly, since £ starts to increase immediately after 
t = 0, we take for initial conditions £(0) = and p(Q) = 1, a choice that leads to C = in equation 
(|22p . As a consequence, Eq. (j22p can be solved for p in the form 



It is clear that £ attains its maximum value £, ma x = 3/2 at r = 3, and returns to its initial value 
£ = at t = 6. Then p changes abruptly from p = — 1 to p = 1 and the motion repeats itself, 




(23) 



This is easily inverted to yield 




(24) 
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as shown in Fig. 2. The period of oscillation is tq = 6 in dimensionless units. These quantitative 
results are not trustworthy and the cusps at the turning points will be smoothed out in the actual 
motion, but the qualitative nature of the dynamics seems to be well represented by our description. 




6 12 18 1 

FIG. 2: Periodic oscillation of the surface of the liquid dielectric near the bottom edge of the capacitor in 
terms of dimensionless quantities in the constant voltage case. 

B. Far from the edge: a quantitative description 

Let us suppose now that Lq ~ L/2, where L is the length of the capacitor, so that the liquid 
starts its motion well inside the capacitor. In this case our appoach is expected to be fairly realistic, 
and furnish quantitative results that can be taken seriously. If the fluid is at rest when the voltage 
is turned on, the initial conditions are x(0) = and x(0) = 0, which are equivalent to £(0) = a 
and £(0) = according to equation (|13p . In this case the constant of integration is given by 
C = a 3 /3 — (1 + a)a 2 /2 and equation ([22]) becomes 

\tV + V® = , V(0 = - ^ - « 2 ) • (25) 

Recalling that p = d^/dr, this equation can be solved in the form 




The appearence of the square root of a cubic polynomial in the above integrand indicates that the 
solution for £(r) will be given by elliptic functions. However, much interesting information about 
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the motion can be obtained without using these not so familiar functions. To this end note that 
equation (|25p is formally equivalent to the one that describes the motion of a particle with zero 
"total energy" in the "potential" V(£). The positivity of the "kinetic energy" £ 2 £ 2 /2 constrains £ 
to the region V(£) < 0. Since x > 0, in terms of £ the physically allowed region is £ > a, as follows 
from Eq.flED. 

The function V(£) has the following properties: (i) V(a) = 0; (ii) V'(a) = —a < which implies 
V(£) < for £ slightly bigger than a; (hi) V(£) — ► oo as £ — > oo; (iv) there is only one positive 
solution to V'(£) = 0, which is readily found to be £ eg = 1 + a and corresponds to x eq = clq. This 
result for £ eg is merely a confirmation of what was found directly from the equation of motion (I14|) . 
Now we have enough information to draw a qualitative graph of the "potential" , which is depicted 
in Fig. 3. The motion is confined to the "potential well" shown in the same figure. 



FIG. 3: Qualitative description of the periodic oscillation of the surface of the liquid dielectric in the constant 
voltage case. The dimensioless variable £ is confined in the "potential well" between £, — a and £ = Cmax- 

The amplitude of the oscillation of the surface of the liquid, that is, the maximum value attained 
by £ is the only positive value £ = £, m ax > a such that V(£) equals zero. In the region £ > £ m ax the 
"potential" is positive and the motion is forbidden. The amplitude ^ ma x is the positive solution to 
the algebraic equation 



V 



a 





(27) 



since £ ^ a. This quadratic equation for £ is solved by 




(28) 
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Typically, as will be seen presently, Lq ^> ao, so that a S> 1. In this case a binomial expansion 
gives the following approximate expression for £, ma x'- 



As a consequence, x max = 2ao- Therefore the surface of the liquid dielectric inside the capacitor 
oscillates periodically between x = and x = x max = 2clq. From (|26|) it follows at once that the 
period of oscillation is 



An expression for T in terms of elliptic integrals is given in the Appendix. 

Unfortunately, under normal laboratory conditions the oscillation of the surface of the liquid 
will be very hard to observe. In order to make ao and to as large as possible, one needs a low- 
density liquid with very high susceptibility. Not only a and b but also their difference should be 
as small as possible, while the voltage should be as high as possible. However, these choices are 
limited, among other things, by the dielectric strength of the liquid. The dielectric strength of 
typical insulators [4] is about 10 7 V/m. The maximum electric field inside the capacitor is easily 
seen to be E max = V/aln(b/a), so that one must have V < 10 7 aln(6/a), with a and b in meters. 
Furthermore, the difference b — a cannot be too small lest capillarity becomes dominant. To give an 
idea of the order of magnitudes involved, consider pure water, for which p ~ 1 g/cm 3 and Xe ~ 79. 
Take a = 1.0 cm, b = 2.0 cm, L = 20 cm, Lq = 10 cm, and V = 5000 volts, a very high voltage. 
Then from (|12p and (|13j) we obtain ao = 0.86 cm, to = 30 IHS cUld Ot — 11.6, SO that ^max — 

13.6. 

It follows that the amplitude of the oscillation is x max ~ 2ao = 1.7 cm. The period of oscillation 
calculated in the Appendix is T = 0.93 s. It should be stressed that these values for the amplitude 
and the period require an extremely high voltage. 




(29) 




(30) 



V. CONSTANT CHARGE CASE 



Making use of © the equation of motion Q becomes 



(31) 
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The same change of variables (|13p . but this time with ao and to replaced by 



_ XeQ 2 Hb/a) I X eQ 2 Hb/a) 

00 4vr 2 e p(6 2 -a 2 )L 2 5 ' V 4ir 2 e p(b 2 - a 2 )L V 1 ' 



reduces equation (|3"Tj) to the dimensionless form 



where a = Lq/clq and the additional dimensionless parameter a is given by 



XeOO 



L 

The same procedure as in the previous case leads to the following first integral for Eq. (|33p : 



(34) 



e P 'e oc, 1 



ln(l + ctO + 



1 



(5 . (35) 



The initial conditions £(0) = a, p(0) = determine C. A qualitative description of the motion can 
be given following the pattern set by the previous case. The shape of the "potential" is exactly the 
same as the one shown in Fig. 3. The motion is periodic and confined to the interval a < £ < £ maa; 
where, now, ^ max is determined by a transcendental equation because of the logarithmic term in 
(|35p . which prevents us from giving an explicit expression for the amplitude of the oscillation. The 
period is also given by a complicated integral involving the square root of a logarithm. 

VI. CONCLUSION 

The motion of a liquid dielectric attracted by a vertical cylindrical capacitor is an unusual 
example of dynamics of a variable-mass system. As we have seen, because the initial mass vanishes 
for motion of the liquid near the lower edge of the capacitor with the most natural initial conditions, 
the velocity change can be qualitatively modelled by a jump determined by the equation of motion 
itself. Assuming that the surface of the liquid remains far from the ends of the capacitor, we exactly 
solved for the motion in the case of constant voltage and found not too complicated expressions 
for the amplitude and the period of the oscillation. If, instead, the charge of the capacitor is held 
constant, the motion is qualitatively the same but quantitative results are different, and demand 
more complicated numerical computations. 
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APPENDIX A: THE PERIOD IN TERMS OF ELLIPTIC INTEGRALS 

According to ([271) the "potential" V(£) has the three real roots £1 = a, £2 = Cmax = 0, £3 = — 7> 
where 



10 1 

a + 3 + 3a\j 1 + 1 ^ 

da ar 



a + 2, 7 = - 



10 1 

3a\/ 1 + 1 K — a — 3 

6a or 



1 + a 



and the approximate equalities hold if a 3> 1. Thus we can write 



(Al) 



v(0 



1 + a 2 



£ 3 - a 3 1 



2 ' 3 

and equation (f30l) for the period takes the form 



(£-«)(/?-£)(£ + 7) 



(A2) 



Jo V — MX Ja 



x dx 



V^V^x) Ja ^/(a;-a)(/3-x)(7 + x) 

We follow [5[ and perform the change of variable x = — u 2 , which leads to 



(A3) 



1 = 2 



(0 - u 2 )du 



20h - 2/ 2 • 



(A4) 



/o v/(/3-q-^)(/3 + 7 -^) 
The additional change of variable u = y/0 — a sin reduces both integrals I\ and I2 to the standard 
form of complete elliptic integrals. Let us first tackle I\: 



h 



'P-a 



du 



v/(/3-q-^)(/3 + 7 -^) V0 + ^ Jo v 7 ! - k 2 sin 2 6 ' 



Tf/2 



de 



that is, 



(A5) 



1 



K(k), k 



1 — a 



where if (&) is the complete elliptic integral of the first kind. Similarly 



(A6) 



13 



h 
so that 



0-a W 2 sin 2 9 d6 



P-a r/ 2 -(l-k 2 sm 2 6)/k 2 + 1/^ 



VP + 7 Jo a/1 - k 2 sin 2 i 



+ 7 -/o 



VT — & 2 sin 2 i 



d0, 



(A7) 



/3 — a 
k 2 ^/W+l 



tt/2 



d0 



tt/2 



1 - k 2 sin 2 dd 



io Vl-k z sm*& JO 
where E{k) is the complete elliptic integral of the second kind. 
After some simplifications, we finally get 



VP + ^r K(k) - E{k) , (A8) 



T = 2^6 1 \ VWTiE{k) 



P + l 



K{k) 



The functions K(k) and E{k) are given by the following power series (0 < \k\ < 1) 



(A9) 



K(k) 



tt/2 



de 



TT 



o vl - k 2 sin 2 i 



x 1+ hr U z + 



1 • 3 
2~ 4 



fc 4 + 



1 • 3 • 5\ 2 
2-4-6 



A; 6 + 



(A10) 



tt/2 



1 - k 2 sin 2 6» d0 



7T 



1 



1 • 3 



l-3-5\ 2 



+ ... 



(All) 



lo 2 I \2J \2 ■ 4/ 3 V2 • 4- 6/ 5 

For the parameters chosen in the main text we have seen that a = 11.6. Then (1A1I) yields 
/3 = 13.6 and 7 = 6.26. With these values the modulus of the complete elliptic integrals is found 
from (|A6P to be k = 0.317. For this k the values for the complete elliptic integrals furnished by 
the above power series are K(k) = 1.6125 and E{k) = 1.5306. Finally, with t$ = 30 ms equation 
01 gives T = 0.93 s. 
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